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LONG-TERM  GOALS 

The  goals  of  this  work  are  to  obtain  a  better  understanding  and  predictive  capability  for  the  processes 
involved  in  deep  ocean  boundary  mixing  caused  by  interactions  of  the  internal  tide  with  sloping 
topography  as  functions  of  key  environmental  parameters,  such  as  ocean  topography,  ambient 
stratification,  latitude,  and  tidal  energy. 

OBJECTIVES 

•  demonstrate  the  principal  mechanisms  that  produce  tidally  induced  boundary  mixing. 

•  quantify  the  magnitude  of  the  mixing  rates  and  parameterize  these  values  into  larger  scale  eddy 
diffusivities  for  momentum  and  density. 

•  calibrate  and  verify  the  model  predictions  and  error  bars  by  testing  against  the  HOME  data  sets 
collected  over  the  Hawaiian  ridge  and  interpret  field  data  in  light  of  results  from  the  numerical 
simulations. 

•  investigate  the  dependency  of  the  mixing  efficiency  as  a  function  of  environmental  conditions 
appropriate  for  other  ocean  sites. 

APPROACH 

The  work  involves  theoretical  development,  numerical  computations,  and  comparison  with  field  data. 
Our  model  formulation  implements  a  non-hydrostatic,  large  eddy  simulation  (LES),  turbulent 
boundary  layer  model  (Slinn  and  Riley,  1998)  on  an/-plane.  Environmental  parameters  are  chosen  to 
represent  conditions  from  the  HOME  site  and  barotropic  tides  are  forced  similar  to  those  present 
during  spring  tides  along  the  Hawaiian  ridge. 
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WORK  COMPLETED 


We  have  completed  a  series  of  numerical  simulations.  Physical  parameters  used  are  taken  from  the 
field  measurements  at  HOME,  e.g.,  N  =  2.2  x  10'3  s’1  (1.25  cph),  the  latitude  is  21.8  N,  giving  a 
Coriolis  parameter  of  2.7  x  10°  s’,  and  the  tidal  period  of  the  M2  tide  is  12.4  hours.  With  our  non- 
dimensionalization,  these  yield  a  background  Richardson  number  of  4.84  and  a  Rossby  number  of  37. 
We  simulate  the  flow  for  5-10  tidal  cycles  to  develop  a  statistical  average  of  mixing  during  different 
phases.  The  x-componcnt  of  the  free  stream  velocity  induced  by  the  baro tropic  tide  is  forced  in  the 
model  to  satisfy  Uoo  =  UMCOs(wt).  Results  from  the  numerical  simulations  are  presented,  followed  by 
an  interpretation  according  to  a  simple  model.  Increased  internal  wave  activity  and  overturning  of 
isopycnals  are  evident  during  phases  of  upslope  flow.  The  signature  of  internal  waves  propagating 
away  from  the  turbulent  region  are  also  evident,  but  decrease  significantly  during  phases  of  down  slope 
flow.  The  velocity  fields  show  deviations  from  the  free  stream  velocity  out  to  heights  of 
approximately  0.35  (70  meters).  These  are  qualitatively  similar  to  time  series  observed  in  the  HOME 
field  data  for  temperature  and  velocity  fields  (not  shown). 

We  are  focusing  on  an  idealized  uniform  slope,  having  horizontal  dimensions  much  larger  than  the 
turbulent  boundary  layer  thickness.  Over  a  tidal  cycle,  flow  external  to  the  boundary  layer  may  travel 
up  and  down  the  slope  on  the  order  of  3  kilometers  (e.g.,  with  UM  ~  20  cm/s  and  a  half  period  of  6.2 
hours).  The  slope  itself  at  the  HOME  Survey  Site  extends  fairly  uniformly  for  approximately  15  km. 
We  are  modeling  a  domain  with  dimensions  of  200  m  on  each  side.  In  the  simulations  presented 
below  the  length  scales  are  chosen  to  be  Lx  =  L  =  Lz  =  200  m  (/.<?.,  non-dimensional  lengths  are  Lx  = 

Ly  =  Lz  =  1.0).  The  typical  turbulent  boundary  layer  thickness  we  observe  is  on  the  order  of  50  m,  and 
the  larger  turbulent  eddies  in  the  domain  are  on  the  order  of  20  m  in  diameter.  In  these  simulations  we 
used  a  grid  resolution  of  129  x  129  x  150  grid  points  (~2.5  million  points).  With  the  grid  clustering 
near  the  boundary  we  resolve  eddies  down  to  approximately  1-2  m  in  diameter  while  smaller  turbulent 
structures  are  modeled  with  the  turbulent  eddy  viscosity  model. 

RESULTS 

Two-dimensional  cross  sections  of  the  non-dimensional  density  field  on  the  x-z  plane  located  at  y0= 
Ly/2  during  three  different  phases  of  the  tide  are  shown  in  the  left  column  of  Figure  1.  The  second 
column  of  panels  in  Figure  1  presents  the  density  field  from  an  “end  view”  on  an  y-z  plane  located  at 
x0=  Lx/2  at  the  same  times.  Together  these  panels  illustrate  the  most  significant  feature  observed  in  the 
simulations:  strong  buoyant  plumes  rise  from  the  boundary  layer  during  a  certain  phase  of  the  tidal 
cycle  that  lead  to  vigorous  turbulent  mixing  in  the  benthic  boundary  layer  for  approximately  0.3  T  (4 
hours).  The  third  column  of  Figure  1  shows  contours  of  the  kinetic  energy  dissipation  rate  nonnalized 
by  the  kinematic  viscosity  and  the  buoyancy  frequency  squared  (si vN2 )  in  the  same  x-z  plane.  The 
right-most  column  of  Figure  1  also  shows  the  dissipation  rate  from  a  “top-view”  in  an  x-y  plane 
parallel  to  the  bottom  boundary  located  at  a  height  zG=  0.05  Lz.  Higher  values  indicate  the  cores  of 
eddies,  the  strong  shear  layer  at  the  boundary,  and  other  coherent  vortical  structures.  The  top  row  of 
panels  is  shown  just  after  a  time  of  maximum  upslope  flow  (i.e.,  t  ~  1)  as  mixing  begins  to  develop 
into  turbulence  and  the  two-dimensional  buoyant  plumes  are  released  from  the  lighter  fluid  near  the 
wall.  The  middle  row  of  panels  is  shown  during  a  phase  of  strong  mixing,  approaching  the  end  of  the 
phase  of  up-slope  flow.  The  bottom  row  of  panels  is  taken  at  a  time  of  strong  stratification,  during  a 
phase  of  down-slope  flow  when  mixing  has  ceased. 
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Figure  1.  Density  fields  (left)  and  dissipation  rate  contours  (right)  during  phases  of flow  instability 
(top  row),  turbulent  mixing  (middle  row)  and  strong  stratification  (bottom). 


The  key  physical  process  uncovered  in  the  simulations  is  that  there  is  a  strong  asymmetry  during  tidal 
phases  caused  by  interactions  between  the  boundary  layer  and  the  background  stratification.  When  the 
flow  is  upslope  heavier  fluid  is  moving  up  from  below.  The  fluid  previously  located  near  the  wall  is 
lighter  than  the  fluid  coming  to  replace  it,  but  the  friction  or  turbulence  near  the  wall  prevents  the 
lighter  fluid  from  leaving  as  quickly  as  the  new  heavier  fluid  is  approaching.  This  leads  to  a  situation 
where  the  lighter,  slower  moving  fluid  is  retained  beneath  the  denser  water  external  to  the  boundary 
layer.  This  unstable  density  structure  greatly  enhances  turbulence  production  near  the  wall.  The 
turbulence  in  turn  produces  a  thicker  boundary  layer  of  slower  moving  intermediate  density  mixed 
fluid.  The  upslope  flow  persists  for  approximately  6  hours  and  turbulence  is  particularly  evident  for 
the  latter  3  hours  when  the  flow  is  decelerating.  The  buoyancy  period,  external  to  the  boundary  layer, 
is  approximately  0.8  hours,  and  hence,  the  stratified  turbulence  responds  in  a  relatively  short  time  scale 
compared  to  the  duration  of  the  flow  conditions,  and  yet  far  from  an  instantaneous  adjustment,  as 
might  be  anticipated  for  buoyant  fluid  with  a  larger  density  difference.  Soon  after  the  onset  of  mixing 
the  intermediate  density  mixed  water  produced  in  the  boundary  layer  is  also  comparatively  light  to  the 
deeper  water  still  flowing  upslope.  The  process  repeats  over  and  over  during  phases  of  upslope  flow 
causing  a  quasi-steady  turbulent  layer  with  weak  or  unstable  stratification  near  the  boundary.  We 
observe  continuous  production  and  release  of  buoyant  plumes  of  water  from  the  boundary  layer  for 
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approximately  3  hours,  during  each  cycle  of  upslope  flow  for  these  conditions.  The  turbulent  layer 
becomes  on  the  order  of  50  -  80  m  thick  in  these  simulations. 

As  the  flow  reverses  downslope,  other  dynamical  processes  become  more  important  and  the  near 
boundary  fluid  containing  less  inertia  responds  to  the  barotropic  pressure  gradient  first  and  moves 
down-slope  while  the  flow  exterior  to  the  boundary  layer  is  still  directed  upslope.  This  final  phase  of 
the  mixing  cycle  leads  to  the  strongest  (shear  induced)  turbulence  close  to  the  boundary.  But  the 
turbulent  layer  soon  collapses,  because  frictional  forces  in  the  boundary  layer  during  the  down-slope 
phase  of  flow  have  the  opposite  effect  from  the  upslope  phase.  Now  buoyancy  forces  act  to  inhibit 
turbulent  production  as  lighter  stratified  fluid  from  above  is  moved  above  the  denser,  slower  moving 
boundary  layer  water  that  is  retarded  by  friction  near  the  wall.  This  eventually  leads  to  much  stronger 
stratification  (by  a  factor  of  2  to  10)  across  the  boundary  layer  than  in  the  quiescent  fluid.  This  stable 
situation  is  sufficient  to  inhibit  turbulent  production  near  the  wall  and  reorganize  the  boundary  layer 
into  a  parallel  shear  flow.  Mixing  ceases  and  remains  negligible  for  approximately  6-8  hours  until  the 
flow  again  returns  upslope.  It  takes  approximately  2  hours  of  upslope  flow,  for  the  specified  physical 
parameters  before  the  turbulence  recommences  in  earnest.  Similar  behavior  is  suggested  by  the  field 
data  from  HOME  and  a  related  asymmetry  between  induced  mixing  on  flood  and  ebb  tides  in  estuaries 
has  been  observed  previously,  e.g.,  Simpson  et  ah,  1990. 


Figure  2.  Illustrative  sketch  of  boundary  layer  production  of  buoyancy  flux.  Left  panel  shows  a 
typical  boundary  layer  profile  that  could  be  idealized  as  a  thinner  layer  of  fluid  moving  at  about  half 
the  velocity  of  the  free  stream.  For  the  tidal  flow  used  here  with  an  upslope  excursion  of  3  km, 
therefore  the  free  stream  would  move  vertically  230  m,  while  the  slower  moving  layer  would  only 
transverse  115  m.  This  would  allow  a  layer  of  lighter  fluid  to  move  115  meters  above  the  heavier 

bottom  retarded  fluid  in  the  boundary  layer. 

We  repeated  the  simulations  with  different  values  of  the  bottom  slope,  Reynolds  number,  and  Prandtl 
number.  The  distance  from  the  boundary  that  the  turbulence  reached  was  less  for  the  other  slopes.  For 
the  9°  slope,  the  mixing  begins  earlier  in  the  cycle  and  persists  for  longer  duration  but  the  plumes  do 
not  extend  as  far  from  the  boundary.  This  is  consistent  with  the  interpretation  (below)  that  there  is 
greater  buoyancy  flux  in  this  case  and  it  is  released  sooner  because  the  vertical  distance  traversed  by 
the  free  stream  velocity  is  twice  as  large  as  for  the  4.4°  case.  No  turbulence  is  produced  in  the  flat 
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bottom  case  indicating  that  the  release  of  tidal  potential  energy,  rather  than  interactions  caused  by  the 
tidal  kinetic  energy  alone  are  responsible  for  the  mixing. 


The  boundary  layer  is  observed  to  become  strongly  restrati  lied  between  mixing  cycles.  The  amount  of 
tidal  potential  energy  that  is  released  to  turbulent  mixing  in  the  benthic  boundary  layer  on  the  slope  is 
directly  related  to  the  boundary  layer  thickness,  or  more  precisely  to  the  momentum  thickness,  0  , 

s  .  f  ..  \ 


defined  as  0{t)  -  J -U— 


0  V 


1- 


U 


dz  .  Here,  u(z)  is  the  horizontally  averaged  instantaneous  upslope 


CO  j 


velocity  profile  and  8  is  the  boundary  layer  thickness,  above  which  the  free  stream  velocity  is  within 
1%  of  Uoo.  For  a  steady  turbulent  boundary  layer  in  an  unstratified  environment,  over  a  flat  plate,  it  is 
common  that  the  momentum  thickness  is  approximately  10%  of  the  boundary  layer  thickness.  A 
laminar  boundary  layer  has  <9  =0.15  5.  It  is  unclear  what  the  relationship  between  9  and  5  in  the 
unsteady  flow  over  a  rough  wall  in  the  benthic  boundary  layer  in  a  stratified  fluid.  Measurements  from 
the  simulations  were  typically  between  7  and  20%  of  8  depending  on  the  phase  of  the  tide. 

The  flux  of  available  potential  energy  produced  by  the  process  of  heavier  upslope  fluid  moving 
faster  on  the  slope  than  lighter  boundary  layer  fluid  is 

0  r\2 


dPEA 

dt 


=  \ypgwdV  =  Axygj  d^zUJt)sm(a)dz  =  Axyg^UK  ^-(f)sin(a)  where  potential  energy  is 
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dp 

dz 


defined  (Winters  et  al,  1995)  as  PE  =  g  [  pzdV .  From  this  we  could  estimate  the  total  mixing 

Jy 

produced  by  this  process  in  the  deep  ocean  over  a  planfonn  area  (Axv)  as 


Mixing  *  Rf/3Axyg 


Q 2  JJ  Q-j 

- M. -J—  a  where  p  is  the  factor  accounting  for  the  angle  between  the  principle 

2  3  dz 

axis  of  the  barotropic  tide  and  the  upslope  direction,  the  small  angle  approximation  is  used  for  sin(j3), 
and  the  factor  UmI  3  accounts  for  the  portion  of  the  wave  period  when  the  tide  is  going  upslope  and 
that  the  average  value  of  Uoo  is  less  than  the  maximum  velocity  Um •  The  greatest  variability  in  the 
formula  will  be  caused  by  the  value  of  0,  which  we  estimate  to  be  in  the  range  of  4  -  40  m  for  different 

dp 

situations,  and  may  be  strongly  dependent  on  the  bottom  roughness,  and  on  — —  which  changes  by 

dz 

approximately  2  orders  of  magnitude  across  the  ocean  depth.  Taking  some  ballpark  numbers  for  these 
quantities,  given  the  large  value  of  Axy,  it  is  simple  to  show  that  the  release  of  tidal  potential  energy  by 
this  process  could  produce  in  the  range  of  0. 1  -  1  TW  of  mixing  in  the  ocean  (cf.  Munk  and  Wunsch, 
1998).  Improved  areal  estimates,  utilizing  methods  of  Egbert  and  Ray  (2001)  are  a  focus  of  our 
ongoing  work. 

IMPACT/APPLICATIONS 

Improved  understanding  of  deep  ocean  mixing  has  potential  benefits  for  society  in  improved  models  of 
global  climate  change  and  turbulence  in  the  ocean  is  relevant  to  ocean  acoustics. 

TRANSITIONS 


The  project  is  in  the  initial  stages.  We  have  developed  the  model  and  are  testing  the  sensitivity  of  the 
oceanic  response  to  environmental  forcing  conditions. 
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The  National  Science  Foundation  is  supporting  a  large  multi-institutional  effort,  led  by  Rob  Pinkel  at 
the  University  of  California,  San  Diego,  Scripps  Institution  of  Oceanography,  on  the  Hawaiian  Ocean 
Mixing  Experiment.  Levine  is  a  participant  in  the  HOME  field  study. 
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